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Abstract. The interpretation of the Meixner-Pollaczek, Meixner and Laguerre polynomials 
as overlap coefficients in the positive discrete series representations of the Lie algebra su(l, 1) 
and the Clebsch-Gordan decomposition leads to generalisations of the convolution identities 
for these polynomials. Using the Racah coefficients convolution identities for continuous 
Hahn, Hahn and Jacobi polynomials are obtained. From the quantised universal enveloping 
algebra for su(l, 1) convolution identities for the Al-Salam and Chihara polynomials and the 
Askey- Wilson polynomials are derived by using the Clebsch-Gordan and Racah coefficients. 
For the quantised universal enveloping algebra for su(2) g-Racah polynomials are interpreted 
as Clebsch-Gordan coefficients, and the linearisation coefficients for a two-parameter family 
of Askey- Wilson polynomials are derived. 



1. Introduction 

The representation theory of Lie algebras and quantum algebras, or quantised uni- 
versal enveloping algebras [6], is intimately linked to special functions of (basic) hyper- 
geometric type, see e.g. [24], [6]. In this paper we consider especially the Lie algebra 
su(l, 1) and its quantum analogue U q (su(l, 1)), and we derive convolution identities for 
certain orthogonal polynomials which occur as overlap coefficients. The idea, which is due 
to Granovskii and Zhedanov [11], see also [23], is to consider (generalised) eigenvectors 
of a suitable element of the Lie algebra which is a recurrence operator in an irreducible 
representation of this Lie algebra. Then there is a relation between these eigenvectors 
and the eigenvectors of this Lie algebra element in the n-fold tensor product of irreducible 
representations of the Lie algebra. From the tensor product decomposition in irreducible 
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representations for n = 2, 3, we obtain identities for these eigenvectors involving Clebsch- 
Gordan and Racah coefficients. In particular, if the overlap coefficients are known in terms 
of special functions, we obtain identities for these special functions in this way. 

For the Lie algebra su(l, I) and the positive discrete series representations a special 
case of this approach is contained in Granovskii and Zhedanov [11], but the result is 
not worked out in detail. Elaborating the method of Granovskii and Zhedanov, Van der 
Jeugt [23] obtains a generalisation of the classical convolution identity for the Laguerre 
polynomials [9, 10.12(41)]. Van der Jeugt [23] also considers the boson Lie algebra b(l), 
a central extension of the oscillator algebra, leading to a generalisation of the convolution 
identity for Hermite polynomials [9, 10.13(38)]. The last identity follows from the previous 
one by a well-known limit transition of Laguerre polynomials to Hermite polynomials, see 
e.g. [15]. 

Apart from the Laguerre and the Hermite polynomials, also the Meixner-Pollaczek, 
Meixner and Charlier polynomials, which all fit into the Askey scheme of hypergeometric 
orthogonal polynomials [3], [15], satisfy a convolution identity of the same form. This is 
a straightforward consequence of the existence of a generating function of a special kind, 
see Al-Salam [1]. It is also known that the Meixner-Pollaczek and the Meixner polynomi- 
als can be interpreted as overlap coefficients in the positive discrete series representations 
of su(l,l), see Masson and Repka [22]. In §3 we show how the method of Granovskii 
and Zhedanov for the two-fold tensor product of positive discrete series representations of 
su(l, 1) leads to a generalisation of the convolution identity for Meixner-Pollaczek poly- 
nomials, from which generalised convolution formulas for Meixner, Laguerre, Charlier and 
Hermite polynomials can be obtained by substitution or by limit transitions. Next using 
the three-fold tensor product representation we obtain a very general convolution iden- 
tity for continuous Hahn polynomials, and similarly for the Hahn and Jacobi polynomials. 
These identities can also be viewed as yielding connection coefficients between two sets of 
orthogonal polynomials in two variables with respect to the same orthogonality measure. 
With this point of view, this result coincides with Dunkl's results [7], [8]. Our derivation 
gives an intrinsic explanation for the occurrence of balanced 4_F 3 -series as connection coeffi- 
cients; they are Racah coefficients. Actually, the interpretation as orthogonal polynomials 
in two variables works in general, and is an intrinsic way to determine the S'-functions in 
[11] and [23] in terms of orthogonal polynomials instead of reducing a triple sum to a single 
sum. 

In §4 we apply the same idea to the quantised universal enveloping algebra L r (? (su(l, 1)) 
and its positive discrete series representations. Due to the non-cocommutativity of the co- 
multiplication, which is needed to define the tensor product representation, the tensor 
product of eigenvectors is no longer an eigenvector in the tensor product representation. 
This can be solved if we restrict to operators related to so-called twisted primitive ele- 
ments in £/q(su(l, 1)), see e.g. [21], [17]. Then the whole machinery works and we obtain a 
generalisation of the Al-Salam and Chihara [2] convolution identity for the Al-Salam and 
Chihara polynomials by considering the Clebsch-Gordan coefficients in the two-fold ten- 
sor product. Going to the three- fold tensor product representations yields a very general 
convolution identity for Askey- Wilson polynomials also involving (/-Racah polynomials, 
and Theorem 4.10 is the key result of this paper. Overlap coefficients are also considered 
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in somewhat more generality in Klimyk and Kachurik [14], but we have to restrict our- 
selves to the twisted primitive elements in order to keep the action in the tensor product 
representations manageable. 

It is interesting to note that in this derivation we have a natural interpretation of 
the continuous Hahn, Hahn and Jacobi polynomials as Clebsch-Gordan coefficients for the 
Lie algebra su(l, 1). Similarly, we have an interpretation of the Askey- Wilson polynomials 
as Clebsch-Gordan coefficients for the quantised universal enveloping algebra U q (su(l 7 1)). 
In §5 we shortly discuss the corresponding result for the quantised universal enveloping 
algebra U q (su(2)), where the g-Racah polynomials then occur as Clebsch-Gordan coeffi- 
cients. This case can be obtained formally from the results for U q (su(l, 1)). Since in the 
dual Hopf *-algebra the so-called zonal spherical elements are known in terms of a two- 
parameter family of Askey- Wilson polynomials, cf. [21], we obtain the explicit linearisation 
coefficients for this subfamily of the Askey- Wilson polynomials. 

It should be remarked that there does not seem to be an appropriate (/-analogue of 
the boson Lie algebra b(l). Either, the Hopf *-algebra structure is lacking, or, as in [13], 
the recurrence in the two-fold tensor product representation seems unmanageable. 

Instead of using generalised eigenvectors we use the spectral theory of Jacobi matrices, 
which we recall briefly in §2. In particular we use this theory to interpret certain recur- 
rence operators in £ 2 (Z + )® n , n = 1,2,3, as multiplication operators in certain weighted 
L 2 -spaces on M n . This approach exploits the theory of orthogonal polynomials, cf. Propo- 
sitions 3.3 and 4.3. 

The notation for (basic) hypergeometric series is the standard one as in Gasper and 
Rahman [10]. Unexplained notions for quantised universal enveloping algebras can be 
found in Chari and Pressley [6]. 

Acknowledgement. We thank Tom Koornwinder for useful comments. The first author 
thanks the Universiteit Gent for its hospitality. 

2. Jacobi matrices and orthogonal polynomials 

We recall some of the results on the spectral theory of Jacobi matrices and the relation 
with orthogonal polynomials. For more information we refer to Berezanskh [4, Ch. VII, 
§1], see also Masson and Repka [22], Klimyk and Kachurik [14]. The operator J acting on 
the standard orthonormal basis {e n \ n G Z + } of £ 2 (Z + ) by 

(2.1) Je n = a n e n+ i + b n e n + a n _i e n _i, a n > 0, b n G K, 

is called a Jacobi matrix. This operator is symmetric, and its deficiency indices are (0, 0) 
or (1, 1). In particular, if the coefficients ci n and b n are bounded, J is a bounded operator 
on £ 2 (Z_|_) and thus self-adjoint. J is an unbounded self-adjoint operator if Xl^Lo a n 1 = 00 
by Carleman's condition. Then eo is a cyclic vector for J, i.e. the span of finite linear 
combinations of the form J p eo, p G Z+, is dense in £ 2 (Z + ). This is the case for all Jacobi 
matrices considered in this paper. 

Assuming this, we can use the same coefficients a n , b n to generate polynomials p n (x) 
of degree n in x by the recurrence relation 



(2.2) xp n (x) = a n p n+1 (x) + b n p n (x) + a„_ip n _i(i), P-i{x) = 0, p Q (x) = 1. 
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By Favard's theorem there exists a positive measure m on the real line such that the 
polynomials p n (x) are orthonormal; 

Pn(x)p m (x) dm(x) = <y n>m . 

The measure is obtained by m{B) = (E(B)e 0l e ) : B Borel set, where E denotes the 
spectral decomposition of the self-adjoint operator J. 

We can represent the operator J as a multiplication operator M x on L 2 (m), where 
M x f(x) = xf(x). For this we define 

A:£ 2 (Z+) ^L 2 (m), (Ae n )(x) = Pn (x), 

then A is a unitary operator, since it maps an orthonormal basis onto an orthonormal basis. 
Note that we use here that the polynomials are dense in L 2 (m), since the self-adjointness 
of J implies that the corresponding moment problem is determined. From (2.1) and (2.2) 
it follows that A o J = M x o A, so that A intertwines the Jacobi matrix J on £ 2 (Z + ) with 
the multiplication operator M x on L 2 (m). 

3. The case su(1, 1) 

The Lie algebra 5u(l, 1) is given by 

[H, B] = 2B, [H, C] = -2C, [B, C] = H. 

There is a *-structure by H* = H and B* = —C. 

The positive discrete series representations Wk of su(l, 1) are unitary representations 
labelled by k > 0. The representation space is i 2 (Z+) equipped with orthonormal basis 
{ e n}{nez+}- The action is given by 

ir k (H) e k n = 2(k + n) e*, 
(3.1) ix k {B) e k n = v / (n + l)(2fc + n) e k n+l , 

n k (C) e k n = -y/nQk + n-Ve*^. 
The tensor product of two positive discrete series representations decomposes as 

oo 

The corresponding intertwining operator can be expressed by means of the Clebsch-Gordan 
coefficients 



p k 2 
c n 2 - 



(3-3) e\= J2 C tnl' h n4\ 

m,n 2 

Later we also use the notation e^^ 2 ^ for ejj to stress the fact that this vector arises 
from the decomposition <S> ^k 2 into irreducible representations. The Clebsch-Gordan 
coefficients are non-zero only if ri\ + ni = n + j, k = k\ + ki + j for j, m, n,2, n G Z + by 
considering the action of H on both sides. We normalise the Clebsch-Gordan coefficients 
by (e fc ,eS 1 ®ef)>0. 

For the above results Vilenkin and Klimyk [24, §8.7] can be consulted. 
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3.1. Clebsch-Gordan coefficients and orthogonal polynomials. The Meixner-Polla- 
czek polynomials are defined by 

(3.4) P(*W) = ^e^ 2Fl (" n ' 2 A A + ";l-e- 2 ^) • 

For A > and < <fi < it these are orthogonal polynomials with respect to a positive 
measure on R, see [15]. The orthonormal Meixner-Pollaczek polynomials 



Pn(x) = P<*W) = J w^o^ W) 



r(n + 2A)' 



satisfy the three-term recurrence relation 



2x sm<pp n (x) = a n p n+ i(x) - 2(n + A) cos(pp n (x) + a n -ip n -i(x), 
a n = y/(n + l)(n + 2X). 

The orthogonality measure for Meixner-Pollaczek polynomials is absolutely continuous. 
Define 

\2A 



w W( X]< f,) = (2sirW>) e(2 ^_ 7r)x | r(A+ . x) ^ 



then 



P 



(x; 4>)p^} (x; 4>) w w (x; 4>)dx = 5 nrn . 



Define the self-adjoint element in su(l, 1); 
(3.5) X (f> = -cos<f)H + B -C. 

Proposition 3.1. A:£ 2 (Z + ) -> L 2 (R, w^ k \x; <f>)dx), e k ^ p£\-; (/>), is a unitary mapping 
intertwining TTf c (X ( j ) ) acting in I 2 with M2 X sin 4> 

on L 2 (R,w (fc) (x; 4>)dx) . 

Here, and elsewhere, M g denotes multiplication by the function g, so M g f(x) = 
g(x)f(x). 

Proof. Use (3.1) and (3.5) to see that ir^X^) is a Jacobi matrix. Next compare the 
coefficients with the three-term recurrence relation for the orthonormal Meixner-Pollaczek 
polynomials to find the result as in §2. □ 

Proposition 3.1 states that v k (x) = Y2n°=oP^ ( x: > ^0 e n * s a generalised eigenvector for 
TTkiX^) for the eigenvalue 2a;sin</>. Next we study the action of in the tensor product 
representation <S> 7Tfc 2 . Recall that A(X^) = 1 <g) X^ + X^ ® 1. 
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Proposition 3.2. T: £ 2 (Z+) ® £ 2 (Z+) -> L 2 (R 2 , w^ kl \ Xl ; (j>)w^\x 2 ; (f>)dx 1 dx 2 ), defined 
by ®e^ 2 \— > Pn^\x\\ (f))p^ 2 2 \x 2 ; <j>) is a unitary mapping intertwining 7T/ Cl ®ixk 2 (A(-X^)) 

Wi/l M 2 (a; 1+a;2 ) s i n ^. 

Proof. This can be seen by using the mapping A of Proposition 3.1 in the second tensor 
factor and solving the resulting three-term recurrence in the first factor. □ 

Proposition 3.2 states that 



v *iM {xuX2)= J2 p^ ) (x 1 ;cf>)p^(x 2 -, ( f>)e k n \®e k 

ni,ri2=0 



2 
ri2 



are generalised eigenvectors for tt^ <S> 7Tk 2 (A(X^)) for the eigenvalue 2{x\ + x 2 ) smcf). 

So T maps the basis e k \ ® e*? onto orthonormal polynomials in two variables. By 
the Clebsch-Gordan decomposition (3.2) there exists another orthonormal basis e k for the 
tensor product representation space. So Ye^ gives another set of orthonormal polynomials 
in two variables in L 2 (R 2 , w^ kl \x\\ (j))w^ k2 \x 2 \ (f>)dxidx 2 ). 

In order to formulate the result we need the continuous Hahn polynomials [15], [16]; 

. n (a + c) n (a + rf) n /-n,n + a + 6 + c + rf-l,a + zx 

3.6 p n (x;a,6,c,d) = i : 3^2 , , , ;1 

n! \ a + c, a + a 

satisfying the orthogonality relations for 3?(a, 6, c, <i) > 



1 

2^ 



T(a + ix)r(6 + ix)T(c — zo;)r(<i — ix) p n (x; a, 6, c, d)p m (x; a, 6, c, cf) dx = 

. T(n + a + c)T(n + a + d)T(n + b + c)T(n + b + d) 



n\(2n + a + b + c + d- l)T(n + a + b + c + d- 1) 



The orthogonality measure is positive for a = c, b = d. 
Proposition 3.3. In L 2 (R 2 , w^ kl \x±; 4>)w^ k2 \x 2 ] (/))dxidx 2 ) we have 

Te k (x 1 ,x 2 ) =p (k \x 1 +x 2 \4>) TeS(xi,x 2 ), 

Teo(xi, x 2 ) = Cpj(x 1 ;ki, k 2 - i{x\ + x 2 ), fci, k 2 + i{x\ + x 2 )), 



C = (-2 sin. 



■ jj\ (2j + 2k t + 2k 2 - l)T(j + 2h + 2k 2 - 1) 



Y(2k 1 +j)Y(2k 2 +j) 



Note that Te^xi, x 2 ) is indeed a polynomial in xi, x 2 . 

Proof. The first statement follows from use of the intertwining of Proposition 3.2 and the 
intertwining of (3.2); 

2(x 1 + x 2 ) smcf)Te k (x 1 , x 2 ) =M 2{xi+X2)sin<p Te k (x u x 2 ) = (T7r fe (X^) e k ) {x u x 2 ), 
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which gives a three-term recurrence relation for Ye^ with respect to n of the same form as 
in Proposition 3.1. Taking into account the initial conditions proves the first statement. 
The prove the second statement we note that for k = k\ + k 2 + j, I = k\ + k 2 + i, 

S ij 8 rnn = {e k ,e l m ) = {Te k ,Te l m ) = [[ pW( Xl + x 2 ; + x 2 ; <f>) 

J Jr 2 

x ^Teo(xi, x 2 )Te l (xi, x 2 ) s jw { - kl \xr, 4>)w^ k2 \x 2 \ (j)) dx\dx 2 , 
by the first statement and Proposition 3.2. Introduce s = x\ + x 2 , t = xi, then we find 

Mm»= / P { n k) (s;<P)p^(s;<P) [ Te^s-t^e^s-t^Xt-^w^Hs-t-^dtds 
Jr Jr 

In case k = I, or % = j, we see that the inner integral must equal the normalised orthog- 
onality measure for the Meixner-Pollaczek polynomials pn(s;cj)), since the corresponding 
moment problem is determined. In case k ^ I, or i ^ j, we conclude that the inner integral 
integrated against any polynomial gives zero, so that it must be zero since the polynomials 
are dense in L 2 (IR 2 , w^ kl \xi; 4>)w^ k2 \x2] 4>)dx\dx2)- So we get 

5 ljW ^(s; 0) = e^-) ^ 23111 ^ 1 ^ 2 Te fc (t, s - t)Te l (t, s - t) 

x r(/ci + it)T{k 2 -is + zt)r(/ci - it)T(k 2 +is- it) dt. 

Apply T to (3.3) for n = to see that Te^ (t, s — t) is a polynomial of degree j in t. Hence, 
Teg(t, s — t) is a multiple of a continuous Hahn polynomial of degree j with the parameters 
as in the proposition. 

The value of the constant follows from comparing the squared norms up to a sign. 
The sign is determined from the condition on the Clebsch-Gordan coefficients. This implies 
< (Te^, TeQ 1 <g> e k2 } and using the first two parts of the proposition and Proposition 3.2 
shows that the sign of C follows from the sign of a double integral of two orthogonal 
polynomials. Only the integral over x 2 is relevant, and the sign of C equals the sign of 
the leading coefficient of the continuous Hahn polynomials viewed as a polynomial in x 2 , 
which is (— 1) J . □ 

So we can now apply T to (3.3) to find, k = k\ + k 2 + j, 
(3-7) £ Cil^p^Hx^p^Hx^) = Cp^(x 1 + x 2 ; ( f>) 

n 1 +n 2 =n+j 

x Pj(xi, hi, k 2 - %{x\ + x 2 ), ki, k 2 + i(x 1 + x 2 )). 

The Clebsch-Gordan coefficients remain to be determined, and this can be done from this 
formula, see [23]. They can be expressed in terms of 3F2-series, which are known as Hahn 
polynomials. Using the Hahn polynomials defined by 

~ , , , TX ^ / — n, n + a + b + 1, — x 
(3.8) Q n (x;a,b,N) = 3 F 2 ( Q+1 _ N ^ 
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for N G Z + , < n < N, we have, with k = k\ + k 2 + j, n\ + n 2 = n + j, 




n\m\n 2 \j\ (2/ci + 2k 2 + 2j) n (2k 2 ) J (2k 1 + 2k 2 + j - 1) 3 

x (n + j)!Qj(ni;2/ci - 1, 2fc 2 - l;n + j), 



(2k 1 ) ni (2k 2 ) n2 (2k 1 ) J 



see [24, §8.7] for another proof. 

Using this in (3.7) gives an identity in a weighted L 2 -space, but since it is a polynomial 
identity it holds for all xi, x 2 . Simplifying proves the following theorem. 

Theorem 3.4. With the notation for continuous Hahn, Meixner-Pollaczek and Hahn poly- 
nomials as in (3.4), (3.6) and (3.8) the following convolution formula holds: 



Remark 3.5. (i) The case j = gives back the convolution identity for the Meixner- 
Pollaczek polynomials, see e.g. [1, §8], [2]. The case n = gives another convolution 
identity for Meixner-Pollaczek polynomials, since the Hahn polynomial reduces to a sum- 
mable 2-Fi-series. 

(ii) Note that the polynomials on both sides of the formula in Theorem 3.4 are orthog- 
onal polynomials in two variables for the space L 2 (R 2 , w^ kl \xi; (f))w( k2 \x 2 ] <j))dxidx 2 ), so 
we have proved a connection coefficient formula for these polynomials. The dual connec- 
tion coefficient formula follows from the orthogonality of the Clebsch-Gordan matrix, or 
equivalently, from the orthogonality relations for the dual Hahn polynomials. 

(iii) Theorem 3.4 shows that the continuous Hahn polynomials have an interpretation 
as Clebsch-Gordan coefficients for su(l, 1). Using the generalised eigenvectors we formally 
have, cf. (3.3), 

v kl,k2 (xi, x 2 ) = ^2 Cpj(xi; ki, k 2 - i(xi + x 2 ), ki, k 2 + i(x\ + x 2 )) v k (xi + x 2 ), 
k 

with C as in Proposition 3.3. The dual relations can be written using the orthogonality 
measure for the continuous Hahn polynomials. 

Recall the definition of the Laguerre polynomials L^\x) = (a + l) n /n! i-Fi(— n; a + 
l;x), the Jacobi polynomials P^ L ' b \x) = (a + l) n /n! 2 Fi(—n, n + a + b + l;a+l;(l—x)/2) 
and the Meixner polynomials M n (x; /?; c) = 2 Fi(— n, —x; (3; 1 — c -1 ). 

Corollary 3.6. (i) ([23]) The Laguerre polynomials satisfy the following convolution iden- 
tity 




Qj(l; 2fci - 1, 2k 2 -l,n + j) P[ kl) On; 0) P^ifa; 4>) = 



Pi kl+k2+i) {x! + x 2 ; 4>) pjixt; h, k 2 - i{x x + x 2 ),h, k 2 + + x 2 )). 



n+j 




E a ^ n + ^ L i a) M L n +J -iM 



X 
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(ii) The Meixner polynomials satisfy the following convolution identity 



( c -i _ ju l ^ n+ 'L Qj(h a-l,b-l,n + j) M,(*i; a; c) M n+J _ z (x 2 ; 6; c) = 

M n (xi + x 2 - j;a + b + 2j; c) (-xi - x 2 ) j Qj(x 1 ;a - 1,6 - l,xi + x 2 ). 



(n + j) 



Proof. The first case follows from the limit transition of the Meixner-Pollaczek polyno- 
mials to the Laguerre polynomials; lim^o P^ a+1 ^ 2 \— 2x/4>; (j)) = L^\x). In this limit 
transition the continuous Hahn polynomials tend to the Jacobi polynomials. 

The second case follows from the substitution = lnc/2z, and replacing x\ and x 2 
by %k\ + %x\ and ik 2 + ix 2 . For this substitution the continuous Hahn polynomials go over 
into the Hahn polynomials. □ 

Remark 3. 7. (i) The case j = in both formulas gives back the convolution identities for 
the Laguerre and Meixner polynomials, see e.g. [1], [2], [9, 10.12(41)], and the case n = 
gives another convolution identity for the Laguerre and Meixner polynomials. Again these 
formulas can be viewed as connection coefficient formulas for orthogonal polynomials in 
two variables. 

(ii) The identities of Corollary 3.6 can be obtained by considering the action of X = 
—H + B — C in the representations nk and ir^ ®^k 2 f° r the Laguerre case, see [23], and by 
considering the action of X c = — (1 + c)/2^/cH + B — C, < c < 1, in the representations 
7Tfc and 7Tfc 1 ® 7Tfc 2 for the Meixner case. The limit case c j 1 in the Meixner result gives 
the Laguerre result. In this case we can interpret the Jacobi and Hahn polynomials as 
Clebsch-Gordan coefficients, cf. Remark 3.5 (iii). 

(iii) Corollary 3.6(h) is equivalent to Theorem 3.4 by the same substitution. Theorem 
3.4 can also be obtained from Corollary 3.6(i) by a double application of the Mellin trans- 
form. For this we have to use that the Laguerre polynomials are mapped onto Meixner- 
Pollaczek polynomials, cf. [18, §3], and that the Jacobi polynomials are mapped onto the 
continuous Hahn polynomials, cf. [16, (3.4) with Y{j3 — iX) replaced by Y{j3 + 

The other hypergeometric orthogonal polynomials satisfying a convolution identity 
are the Charlier and Hermite polynomials, cf. [1], [2]. These identities can be obtained 
by taking the appropriate limits from the Meixner polynomials to the Charlier polynomi- 
als and from the Laguerre polynomials to the Hermite polynomials, cf. e.g. [15]. The 
Hahn polynomials tend to Krawtchouk polynomials and the Jacobi polynomials tend 
to Hermite polynomials. We use the notation K n (x;p, N) = 2 Fi(—n,—x;—N;p~ 1 ) for 
Krawtchouk polynomials, C n (x;a) = 2 Fo(—n,—x;—;a~ 1 ) for Charlier polynomials and 
H n (x) = (2x) n 2 F (—n/2, —(n — l)/2; — ; — x~ 2 ) for the Hermite polynomials. 



Corollary 3.8. (i) ([23]) The Hermite polynomials satisfy the following convolution iden- 
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tity 



n +j „2 „l frn+j-l 



E K i ft ^2 ^ + j)^H l (x) {n U + . _ , H n+j _ t (y) 



(a 2 + 6 2 )(^')/ 2 (b\ j „ fax + by\ TJ ( ay-bx 



(n + j)\ \aj n Wa 2 + b 2 J ° \Va 2 +b 2 J ' 
(ii) The Charlier polynomials satisfy the following convolution identity 

n+j , .x 

E (" 7 J J (/; n + j)3 (x; a) C B+i _, (y; 0) = 

(-l) J (a + (3) n C n {x + y-j;a + f3) (-x - y) J K J (x; x + y). 

Remark 3.9. (i) Again the case j = gives known convolution formulas, cf. [1, §8], [2], 
[9, 10.13(40)]. Corollary 3.8(ii) is derived in a different way in Vilenkin and Klimyk [24, 
§8.6.5]. 

(ii) This time the identities have a similar interpretation, but now we have to use the 
Lie algebra b(l), a central extension of the oscillator algebra, cf. [23]. In particular we can 
now interpret the Hermite and Charlier polynomials as Clebsch-Gordan coefficients. 

3.2. Racah coefficients and orthogonal polynomials. In the tensor product of three 
positive discrete series representations 7T/ Cl <g> %k 2 ® 7Tfc 3 of su(l, 1) we consider the following 
orthogonal bases; 



(o Q\ J(kik 2 )k 12 k 3 )k _ \^ t~ik 12 ,k 3 ,k (k 1 k 2 )k 12 ~ k 3 

v d - y J e n — 2.^1 °ni 2 ,n 3 ,n S 12 <® e n 3 

ni 2 ,n 3 

/Qin^l — \ r~<k 1 ,k 2 ,k 12 y~ik 12 ,k 3 ,k ki ~ k 2 ft k 3 

W-IV) - 2^ °«i,«2,«l2%,n3,n ^i^^^^s' 



ni ,n 2 ,n 3 ,n 12 



and 

/Oil"! Jki(k 2 k 3 )k 23 )k _ r-ik 1 ,k 23 ,k ii « (fefejtj 



(3-12) = V ni n 2 2 3 3 Cni n ^ ® ej* ® e* 3 . . 



ni,n 2 ,n 3 ,n 23 



Here we use the extended notation ei^ lk2 ^ k for the basis of the tensor product decomposition 
to keep track of how the decomposition is obtained. 

These bases are connected by the Racah coefficients, which leads to an intertwiner for 
the action of su(l, 1). The Racah coefficients are defined by 

/q i o\ „((fcifc 2 )fci2A;3)A; _ jjkx,k 2 ,k\ 2 Ak\{k 2 k 3 )k 23 )k 

W- 1 - 6 ) e n - 2^ U k 3 ,k,k 2 3 e n 
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In the previous formulas the following constraints hold 

ki2 = h + k 2 + ji2, k 23 = k 2 + k 3 + j 23 , 
k = k 12 + k 3 + j = k x + k 23 + j', j 12 ,j, 323, f e Z+, and j 12 +j= j 23 + j' . 

Thus all above sums are finite sums. 

Recall that (1 <g> A)(A(X (j) )) = l®l®X <t> + l®X < i > ®l+X < i ) ®l®l. The following 
proposition is proved as Proposition 3.2. 

Proposition 3.10. Define the unitary mapping 

Q:l 2 (Z + ) ®l 2 (Z + ) ®1 2 {Z + ) -> L 2 (M 3 ,«; (fcl )(xi;(/))«; (/£2) (x2;(/))w (fc3) (x3;(/))^i^ 2 rfx3) 
by 

6: e* l x ® e* 2 ® e*J ~ p (*i)(x i; </>)^ 2 2) (*2; ^pft^a*; 0), 

i/ien 6 intertwines n kl <S> ^k 2 <8> 7Tfc 3 ((l ® A)(A(X^)) wrat/i M 2 ( Xl+X2+X3 ) sin< ^. 

Remark 3.11. Let A^ fc ^ = A be the unitary mapping defined in Proposition 3.1 and 
r ^(k 1 k 2 ) = x be the unitary mapping defined in Proposition 3.2. Using the identifications 

L 2 (R 3 , w {kl) On; (j))w {k2) (x 2 ; (fyw^ (x 3 ; 0)efoiete 2 efe 3 ) 

= L 2 (R, (fydxx) <g> L 2 (R 2 , «; (fc2 )(;r 2 ; (P)w {ka) (x 3 ; (j))dx 2 dx 3 ) 

= L 2 (R 2 , <p)w {k2) (x 2 ; 0))efoitte 2 ) <g> L 2 (R, «;^ 3) (x 3 ; 0)efo 3 ), 

we have = A^ 1 ) <g> T^ 3 ) = T( fclfc2 ) <g> A^. Hence, for the orthogonal bases on the 
right hand side of (3.9) and (3.12) we have 

9 e (fcifc2)fci2 ^ e *s = ^^e^ 2 ^ 12 ) (A (fc3) e^ 3 3 ), 
©e^ ® e^f 3)fc23 = (A (fcl) e^) (r^ ka) e^f 3 ^. 

And the right hand sides are known from Propositions 3.1 and 3.2 in terms of Meixner- 
Pollaczek polynomials times continuous Hahn polynomials. 

Proposition 3.12. (i) The following expressions hold; 

Q( e (i^k 2 )k 12k3 )k ){xu ^ Xs) = p (k) {xi + X2 + Xs . 0) e(ef lfc2)fcl2fc3)fc )(x 1 , x 2 , x 3 ), 
e(e« fclfe2)fel2fe3)fe )(a; 1 , x 2 , x 3 ) = T^e? 1 ^ 12 ^, x 2 ) T^ k ^ 12ka)k ( Xl + x 2 , x 3 ). 

(ii) TTie following expressions hold; 

e(e^ k > k ^ k )( Xl ,x 2 , x 3 ) = pl k \ Xl + x 2 + x 3 - </>) e(e( fcl(fc2fc3)fc23)fc )(x 1 , x 2 , x 3 ), 
e( e ( fc i( fc2fc 3)fe23)fc )(a , i5a , 2)a , 3) = T ( k ^e {k2ks)k23 (x 2 ,x 3 ) T^ k ^e {klk23)k (x 1 ,x 2 + x 3 ). 
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Proof. Statement (ii) is proved analogously as statement (i). The first statement of (i) 
follows from Proposition 3.10 and the decomposition of the three-fold tensor product, cf. 
Proposition 3.3. 

For the second statement we use (3.9), Remark 3.11 and Propositions 3.3 and 3.1 to 

find 

&(e« klk > )kl > k3)k )(x 1 ,x 2 ,x 3 ) = 

ni2+n 3 =j 

The sum can be evaluated as (r^^ef 12fes)fe ) (x 1 + x 2 ,x 3 ) by (3.7). □ 

Next we apply to (3.13), then it follows from Proposition 3.12 that we can divide 
both sides by the Meixner-Pollaczek polynomial of degree n. Since O is unitary we obtain 
the Wigner-Eckart theorem, stating that the Racah coefficients in (3.13) are independent 
of n. So we can restrict to the case n = of (3.13) before applying 6 without loss of 
generality. We obtain 

(3.15) U tlf£l 2 (x^^ef 2 ^ 23 ) (x 2 , x 3 ) (T^ k ^e { hlk23)k ) (x u x 2 + x 3 ) = 

The Racah coefficients remain to be determined, and this can actually be done from 
(3.15), see [23]. One can either copy the expression [23, (4.8)], or use the limit q f 1 of the 
expression for the (/-Racah coefficient given in Proposition 4.9. Both lead to the following 
expression of the Racah coefficients in terms of balanced 4i 7 3-series; 



323 



k,MM2 _ fj+Ji2\ (2k 2 ) J12 (2k 3 ) J (2k 1 + 2k 2 + 2k 3 +j+j 12 - V 



(3-16) U kMa ^ j 



323 



(2k 3 , 2k 2 + 2k 3 +j 23 - l) j23 (2k 2 + 2k 3 + 2j 23 ] 



j'\(2k u 2k 23 ,2k 1 + 2k 23 +f - j 23 \(2k 2 ,2k 3 ,2k 2 + 2k 3 + j 23 - 1) J23 ^ 1/2 



jl(2k 12 , 2k 3 , 2k 12 + 2k 3 +j- I) 3 j 12 \(2k u 2k 2 , 2h + 2k 2 + j 12 - % 

2ki + 2k 2 + ji2 - 1, 2k 2 + 2k 3 + j 23 - 1, -j 12 , -j 23 



x 



4 3 V 2k 2 ,2h + 2k 2 + 2k 3 +j+j 12 -l,-j-j 12 

with the convention (3.14). 

The Racah coefficients can be rewritten in terms of the Racah polynomials defined by 

/o-,^ t-> / \ / \ n f-n,n + a + /3 + l,-x,x + <y + 5 + l 

(3.17) R n {\(x);a,[3,~f,5) = 4 F 3 [ ;1 

\ a + 1, p + + 1, 7+1 

where X(x) = x{x + 7 + 5 + 1), one of the lower parameters equals — N, N e Z + and 
< n < N, cf. [15]. The orthogonality relations for the Racah polynomials follow from 
the fact that the Racah coefficients form a unitary matrix. 

So we obtain the following theorem by simplifying (3.15) using s = x\ + x 2 + x 3 and 
the explicit expression (3.16). 
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Theorem 3.13. The continuous Hahn polynomials satisfy the following convolution iden 
tity 

n+j 

E 



1=0 



j + n\ (2k 2 ) n {2k 3 ) j {2k 1 + 2fe 2 + 2k 3 + j + n - 1) 



n J (2fe 3 ) z (2fe 2 + 2fe 3 + I - l)i{2k 2 + 2k 3 + 2/) i+n _, 
x Ri(X(n); 2k 2 - 1, 2fe 3 -1,-j-n- 1, 2fei + 2fe 2 + j + n - 1) 
x p n+J _/(xi; fei, fe 2 + fe 3 + J - is, fei, fe 2 + fe 3 + / + is) 
x pz(x 2 ; fe 2 , fe 3 - i(s - ffi), fe 2 , fe 3 + i(s - xt)) 
= p n (xi; fei, fe 2 - + x 2 ), fei, fe 2 + + x 2 )) 

x pj(xi + x 2 ; ki + k 2 + n, fe 3 - is, ki + k 2 + n, fe 3 + is), 

with the notation as in (3.6), (3.17). 

Remark 3.14- (i) Theorem 3.13 can be considered as a connection coefficient formula be- 
tween two systems of orthogonal polynomials for the orthogonality measure 

T(ki + ixi)T(k 1 -ix 1 )T(k 2 +ix 2 )T(k 2 -ix 2 )T(ks+i(s-x 1 -x 2 ))T(ks-i(s-xi-x 2 )) dxidx 2 

on R 2 . This follows from substituting s = x\ + x 2 + x 3 in the weighted L 2 -space of 
Proposition 3.10 and leaving out the integration with respect to s, which can be done by 
Proposition 3.12 and the Wigner-Eckart theorem. 

(ii) Theorem 3.4 can be obtained as a limit case of Theorem 3.13 by letting fe 3 — > oo 
and using the limit transition of the continuous Hahn polynomials to the Meixner-Pollaczek 
polynomials, see e.g. [15]. Note that Theorem 3.4 is used in the derivation of Theorem 
3.13. 

(iii) Application of © to (3.9)- (3.12) gives results which are immediately derivable 
from Theorem 3.4. 

Corollary 3.15. (i) ([23]) The Jacobi polynomials satisfy the convolution identity 
( j + n\ (b + l) n (c+ l)j(a + b + c + j + n + 2) t 



1=0 



n J (c+l)i(b + c + l + l)i(b + c + 2l + 2) j+n _i 
x Ri(X(n);b, c, -j - n - 1, a + b + j + n + 1) 



x j££r*«>(i - 2 Xl) (i - xi) >pf* ( 1 7^- 2 ^ ) 

= (n + i 2 )"P<-»> ( ^Z-\) p(«+»+=»+i,c) (1 _ 2{xi + X2)) 

\Xi +- X 2 J 



(ii) The Hahn polynomials satisfy the following convolution identity 

n+j 

E 



" +J 'j + n\ (a+l) n+j _i(b+l) l (b + l) n (c+l) j (a + b + c + j + n + 2) l 



1=0 



I J (a + l) n (c+l)i(b + c + l + l)i(b + c + 2l + 2) j+n _i(a + b + 2n + 2) j 
x Ri(X(n); b, c, -j - n - 1, a + b + j + n + 1) 

x (/ - s) n+j -i Q n +j-i(x\;a, b + c + 21 + 1, s - I) {x\ - s)i Qi(x 2 ; b,c,s- x\) 
(—xi - x 2 ) n Q n (xi, a, b, x\ + x 2 ) (n - s)j Qj{x\ + x 2 - n; a + b + 2n + 1, c, s - n) 
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with the notation (3.8), (3.17). 

Proof. The first result follows from the limit transition of the continuous Hahn polynomials 
to the Jacobi polynomials. Replace Xi by sxi and let s — > oo. The second result follows 
by a similar substitution as in the proof of Corollary 3.6(h). □ 

Remark 3.16. (i) Similar as in Remark 3.7(iii) we have that Corollary 3.15(h) and Theorem 
3.13 can be obtained from each other by formal substitution. Theorem 3.13 can be obtained 
from Corollary 3.15(i) by a double application of the Mellin transform. Moreover, Corollary 
3.15 can be proved as Theorem 3.13 by analysing the action of X and X c , cf. Remark 
3.7(h), in the three-fold tensor product. 

(ii) Dunkl [7, Thm. 4.2, Prop. 5.4], [8, Thm. 1.7] has obtained Corollary 3.15, and 
hence Theorem 3.13, by a different method. Dunkl [7] obtains the two- variable Hahn 
polynomials by judiciously guessing solutions for a certain difference equation arising from 
the representation theory of the symmetric group. By symmetry considerations there are 
more solutions of this type, and the connection coefficients can be calculated in terms of 
balanced 4_F3-series. The derivation in this paper gives an intrinsic explanation for the 
occurrence of the Racah polynomials as connection coefficients. See also Dunkl [7] , [8] for 
the orthogonality relations for these two- variable Hahn and Jacobi polynomials for suitable 
restrictions on the parameters. 

We do not obtain extensions of Corollary 3.8 in this way. For /ci,/c2,&3 — > oo in 
Theorem 3.13 we obtain the same result. This is also explained by the fact that the Racah 
coefficients for the Lie algebra b(l) are of the same form as the Clebsch-Gordan coefficients, 
cf. [23]. 

4. The case U q (su(l,l)) 

Let U q (sl(2,C)) be the complex unital associative algebra generated by A, B, C, D 
subject to the relations 

A 2 - D 2 

(4.1) AD = 1 = DA, AB = qBA, AC = q~ 1 CA, BC-CB= -. 

q — q~ 

It is a Hopf algebra. We are only concerned with the comultiplication, which is defined by 

A(A)=A®A, A(B) = A&B + B® D, 

(4.2) 

v ; A(C) = A®C + C® D, A{D)=D®D 

on the level of generators and extended as an algebra homomorphism. There are several 
possible *-structures on U q {s\.{2, C)), and we take 

A* = A, B* = -C, C* = -B, D* = D, 

and the corresponding Hopf *-algebra is denoted by U q (su(l, 1)). 



CONVOLUTIONS FOR ORTHOGONAL POLYNOMIALS 



15 



The positive discrete series representations n k of U q (su(l, 1)) are unitary represen- 
tations labelled by k > 0. They act in £ 2 (Z + ) and the action of the generators is given 
by 



n k (A) e k n = q k+n e k n: 

n k (C)e k n = q ^- k - n ^- J >-<■>;, ,, 



(4.3) n k y^,e n -qi — 



Kk{B)e n = q _^ e n+1 . 

Note that 7r k (D) is an unbounded operator, but that ir k (A), 7r k (B), ir k (C) £ B(£ 2 (Z + )). 
The operators that we consider are bounded. 

Recall that the tensor product of two representations are defined by use of the comul- 
tiplication. The tensor product of two positive discrete series representation decomposes 
as for the Lie algebra su(l, 1); 

oo 

(4-4) n kl <g> n k2 ^ ^k 1 +k 2 +j 

So there exists a unitary matrix mapping the orthogonal basis e k \ <E> e k 2 onto e kl+k2+ i 
intertwining the action of U q (su(l, 1)). The matrix elements of this unitary mapping are 
the Clebsch-Gordan coefficients; 

oo 

(4-5) «£= £ C k n ^ n e k n \®e% 

ni,ri2=0 

where k = k\ + k2 +j for j e Z + . The sum is finite; n\-\-ni = n + j. The Clebsch-Gordan 
coefficients are normalised by (e§, e^ 1 <g> e^ 2 ) > 0. 

These results can be found in Burban and Klimyk [5] and Kalnins, Manocha and Miller 
[13]. See Chari and Pressley [6] for general information on quantised universal enveloping 
algebras. 

4.1. Clebsch-Gordan coefficients and orthogonal polynomials. For this section we 
need the Askey- Wilson polynomials and the Al-Salam and Chihara polynomials, which are 
a subclass of the Askey- Wilson polynomials. The Askey- Wilson polynomial is defined by 

(4.6) p m (cos6>;a, b, c, d\q) = a (ab, ac, ad; q) m 4^3 I ^ ac ad iQ,Q) 

and it is symmetric in its parameters a, 6, c and d, see [3]. The Al-Salam and Chihara 
polynomials are obtained by taking c = d = in the Askey- Wilson polynomials; 

/ -m ae i9 ae ~i9 \ 

(4.7) s m (cos(9;a,%) = p m (cos(9;a, 6,0, 0|g) = a m (a6; g) m 3 y? 2 ( a ^ q 5 9» 9 ) - 
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By dm(-;a,b,c,d\q) we denote the normalised orthogonality measure for the Askey- 
Wilson polynomials, which is absolutely continuous on [—1, 1] and has at most a finite 
number of discrete mass points outside [—1, 1]. We put dm(-; a, b\q) = dm(-; a, 6, 0, 0|g) for 
the normalised orthogonality measure for the Al-Salam and Chihara polynomials. Explic- 
itly, let 

(az, a/z, bz, b/z, cz, c/z, dz, d/z; g)^ 

and we use w(z) = w(z;a,b,c,d\q) to stress the dependence on the parameters when 
needed. Let a, 6, c and d be real, or, if complex, appearing in conjugate pairs, and let all 
the pairwise products of a, 6, c and d not be greater or equal than 1. Then the Askey- 
Wilson polynomials p n (x) = p n (x; a, 6, c, d\q) satisfy the orthogonality relations 

/ Pn(cOs9)p m (cOs9)w(e ld )d9 + p n (x k )p m (Xk)w k = Sn,rnh n , 

2nh Jo h 

^ g-j ^ (1 — q n ~ 1 abcd) (q,ab,ac,ad,bc,bd,cd;q) n 

(1 — q 2n ~ 1 abcd) (abcd;q) n 
(abed; q)^ 



h 



(q, ab, ac, ad, be, bd, cd; q) c 



The points x^ are of the form \{eq k + e~ 1 q~ k ) for e any of the parameters a, b, c or d 
with absolute value greater than 1; the sum is over k G Z + such that \eq k \ > 1 and Wk is 
the residue of z \— > w(z) at z = eq k minus the residue at z = e~ 1 q~ k . So the normalised 
orthogonality measure dm(-; a, b, c, d\q) can be read off from (4.8), see Askey and Wilson 
[3] or [10]. 

Let S m (x; a, b\q) = s m (x; a, b\q) / a/ (q, ab; q) m denote the orthonormal Al-Salam and 
Chihara polynomials, which satisfy the three-term recurrence relation 

2x S n (x) = a n+1 S n+1 (x) + q n (a + b) S n (x) + a n S n -!(x), 

(4.9) . 

a n = y/il-abqn-^il-q"). 



We now define 



(4.10) Y s = q x l 2 B - q-^C + ^^( A ~ D ) e U q (su(l, 1)). 

Then Y S A is a self-adjoint element in U q (su(l, 1)) for s G R\{0}, or s G T. Y s is twisted 
primitive element, i.e. A(Y S ) = A <g> Y s + Y s <g> D, meaning that Y s is very much like a Lie 
algebra element. 

We also use the notation p(x) = (x + x~ l )/2 = p(x~ l ) for x ^ in this section. 

Proposition 4.1. Let A:£ 2 (Z + ) — > L 2 (M, dm(-; q 2k s, q 2k / s\q 2 )) be the unitary mapping 
defined by A: e k i— > S n (-; q 2k s, q 2k / s\q 2 ), then A intertwines 7Tk(Y s A) acting in £ 2 (Z + ) with 
2(M x -p(s))/(q- 1 -q). 

Proof. The bounded self-adjoint operator 7Tk(Y s A) is a Jacobi matrix by (4.10) and (4.3), 
and the result follows upon comparing with the three-term recurrence (4.9) for the Al- 
Salam and Chihara polynomials as in §2. □ 
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Proposition 4.1 says that v k (x) = J2^ =0 S n (fj>(x); q 2k s 7 q 2k / s\q 2 ) e k is a generalised 
eigenvector of the self-adjoint operator nk(Y s A) for the eigenvalue 

_ x + x~ x - s - s -1 _ n{x) - n(s) 

*x — Zi — Zi • 

q - q q ~q 

Due to the fact that the co multiplication on U q (su(l, 1)) is less simple than for the 
Lie algebra su(l, 1) it takes a little more effort to determine the action of Y S A in 7r/ Cl ®7Tfc 2 . 
The result can still be phrased using orthogonal polynomials in two variables. 

Proposition 4.2. Define T: £ 2 (Z+) <g) £ 2 (%+) -> L 2 (R 2 , dm{x u x 2 )), where 

dm(xi, x 2 ) = dm(xi; q 2hl w 2 , q 2kl /w 2 \q 2 ) dm(x 2 ; q 2k2 s, q 2k2 /s\q 2 ), x 2 = ^(w 2 ), 

by 

T:e k n \ ®e\\ » S ni (x i; q 2k ^w 2 , q 2k > /w 2 \q 2 ) S n2 (x 2 ;q 2k2 s,q 2k2 / s\q 2 ) 

then T is a unitary mapping intertwining 7T/ Cl <E> ^k 2 (A(Y S A)) with 2(M Xl — /j,(s))/(q~ 1 —q) 
in L 2 (dm(xi, x 2 )). 

Note that T(e*£ <S> e k2 2 ) forms a set of orthogonal polynomials in two variables x\ and 
x 2 for L 2 (R 2 , dm(xi, x 2 )), since the Al-Salam and Chihara polynomial is symmetric in its 
parameters. 

Proposition 4.2 states that the vector 

oo 

w(x 1 -x 2 ) = S ni (/j(x 1 );q 2kl x 2 ,q 2kl /x 2 \q 2 )e kl i ®v k2 (x 2 ) 
m=0 

oo 

= S n M^q 2kl ^2,q 2kl /x 2 \q 2 )S n M^q 2k2 s,q 2k2 /s\q 2 )e k \ ® e k n \ 

n 1 ,n 2 =0 

is a generalised eigenvector of 7T/ Cl ®TVk 2 (A(Y s A)) for the eigenvalue A^. This last observa- 
tion is essentially the way to obtain Proposition 4.2, since A(Y S A) = A 2 <g> Y S A + Y S A <g> 1 
acts as a three-term recurrence operator in e^ 1 <S> v 2 (x 2 ). 

Proof. We use A(Y S A) = A 2 <g> Y" S A + y s A ® 1 and Proposition 4.1 to define for fixed x 2 
the map A : £ 2 (Z+) ® £ 2 -> £ 2 (Z + ) by 

Ao: e£ ® e* 2 2 ^ S n2 (x 2 ; </ 2fc2 s, (j 2 * 2 /*!? 2 ) 
to obtain the recurrence in ni 

Ao((? _1 -q)(^ kl ® 7r k2 (A(Y s A)) + s + s' 1 ^) e k \ ® e k n \ = 

S n2 (x 2 ; q 2k2 s, q 2k2 /s\q 2 ) (q 2 ^ ((s + s^q 2 ^ + \ W2 q 2k > (q' 1 - q)) e k \ 



+ v /(l- g 2n 1+ 2 )(1 _ ? 4 fel+ 2n )e fei +i + ^ (1 _ g2„ l)(1 _ ^+2^-1) e Ji ^ 

Use the explicit expression for A W2 and the three-term recurrence relation (4.9) to obtain 
the result. □ 



We now calculate the action of T e„, which yields another set of orthonormal polyno- 
mials for L 2 (R 2 , dm(xi, x 2 )). 
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Proposition 4.3. Let k = k± + k 2 + j for j G Z + , and x\ = ^(wi) then 

(T ej) (xi, x 2 ) = S n ( Xl ;q 2k s, q 2k / s\q 2 ) (T e k ) (x u x 2 ), 

(T e k ) (xi, x 2 ) = C Pi (x 2 ; ? 2fcl «;i, g^M, ^Ak'), 

C" 1 = ((Cj(/Ci, k 2 ))~ 1 = ^J(q 2 , g4*i, g4*2 } g4fci+4fc 2 +2j-2. g 2) 

Proof. By Proposition 4.2, (4.4) and 

2 ^przy Te n(^i' *a) = T(7r fc (y s A) e*)^!, x 2 ) 

we obtain the three-term recurrence relation as in Proposition 4.1, but with different initial 
conditions. Hence, the first statement follows. 

Since T is unitary we have the orthogonality relations 5 nm 5ki = (Te„, Ye l m ) = 

S n (xi;q 2k s, q 2k /s\q 2 )S m ( Xl ; q 2l s, q 2l /s\q 2 ) / TeJ(xi, x 2 )Tej,(xi, x 2 ) dm(xi, x 2 ), 



by our first observation. As in the proof of Proposition 3.3 we conclude Te^xi, x 2 ) = 
Pj(x2) is polynomial of degree j, k = k\ + k 2 + j, in x 2 satisfying the orthogonality 
relations 

/ Pj(x 2 )pi(x 2 ) dm( Xl , x 2 ) = 5ijdm( Xl ;q 2k s, q 2k / s\q 2 ) 

J x 2 



as measures with respect to functions in the variable x\. 

We now assume for ease of presentation that dm{x\,x<i) is absolutely continuous. 
The general case can be proved similarly, or it can be obtained by analytic continuation 
with respect to s. The measure is absolutely continuous for q 2k2 < \s\ < q~ 2k2 , since 
ki, k 2 > 0. Put x\ = cos6>, x 2 = cosip, then we obtain the explicit expression (4.8) for the 
orthogonality measure; 

13 (q 2 , q^ , q 4k2 ;q 2 )oo ( g 2fc 1 +2fc 2 +2i se ±ifl j q 2k 1+ 2k 2 +2j e ±i6 / s . q 2^ 

for almost all 9. The ±-signs means that we take all possible combinations in the infinite q- 
shifted factorials. Cancelling the (e ±2l6> ; q 2 )oo on both sides and comparing the result with 
(4.8) we see that pj is a multiple of pj(-; q 2kl e l °, q 2kl e~ l6 , q 2k2 s, q 2k2 /s\q 2 ). The constant in 
front follows up to a sign by comparing the squared norms. As in the proof of Proposition 
3.3 the sign of C follows from the normalisation of the Clebsch-Gordan coefficients, and 
now we obtain C > 0. □ 

So we obtain a second set of orthonormal polynomials for L 2 (R 2 , dm(xi,x 2 )) in terms 
of Al-Salam and Chihara polynomials and Askey- Wilson polynomials. 



CONVOLUTIONS FOR ORTHOGONAL POLYNOMIALS 



19 



The convolution formula for the Al-Salam and Chihara polynomials is obtained by 
applying T to (4.5) using the results of Propositions 4.2 and 4.3. The results holds as an 
identity in a weighted L 2 -space, but since it is a polynomial identity it holds for all x±, x 2 ; 
with x\ = fi(wi), X2 = ^(w 2 ), and k = k\ + k 2 + j 

(4.11) C n\ttn S m q 2kl /w 2 \q 2 ) S n2 (x 2 ; q 2k2 s, q 2k2 /s\q 2 ) = 

n 1 +n 2 =n+j 

S n (xi;q 2h s, q 2k /s\q 2 )pj(x 2 ; q 2kl w u q 2kl /w!, q 2k2 s, q 2k2 /s\q 2 ) 

y/(q 2 , g 4fc 1? ? 4fc 2j g4fc 1+ 4fc 2 +2j-2. q 2). 

We have not yet calculated the Clebsch-Gordan coefficients explicitly, but we can now 
use (4.11) to determine C k \'^^ by specialising to a generating function for the Clebsch- 
Gordan coefficients. The result is phrased in terms of g-Hahn polynomials, which are 
defined as follows; 

(Q~ n : Q~ x : cthq n ~^~^ 
aq, q N 

See e.g. [13] for other derivations of the following lemma. 
Lemma 4.4. With n\ + n 2 = n + j we get 

Ctnli 1+k2+J = C Qj(q- 2ni ;q' k1 - 2 , q' k2 ~ 2 , n+j; q 2 ), 
with the constant C given by 

q 2k l{ n-n l){q 2. q *) n+jy /tfki ; g2) ^ (g 4/c 2 j ^ ^ .^2j~ 

yV; Q 2 )n(q 2 ; q 2 ) ni (q 2 ; q 2 )n 2 (q 2 ; q 2 )j(q^+^+^ ; q 2 ) n (q^; q 2 ) j (q^+ 4k ^- 2 ; q 2 )j ' 

Proof. Observe that C^J'n^'n * s independent of s, x\ = fi(wi) and x 2 = n(w 2 ). Specialise 
w 2 = q 2k2 s and w\ = q 2kl /w 2 = q 2kl ~ 2k2 /s, then the Al-Salam and Chihara polynomials 
in the summand on the left hand side of (4.11) can be evaluated explicitly, since the 3<p 2 - 
series reduces to 1. For this choice the Askey- Wilson polynomial on the right hand side can 
also be evaluated explicitly, and we obtain the generating function for the Clebsch-Gordan 
coefficients 



Efik\,k 2 ,k 2n\(k 2 — k\) — 2n 2 k 2 „ni- n 2 V(Q 4kl ^ 2 )nAq 4k2 ;q 2 ) n 2 
ni+n2=n+j VW »9 )m\Q ,q )n 2 



q -2jk 2 -2n(k 1 +k 2 +j) s n-j^ q 4k 1 ^ q 4k 2 ^ ^4fc 2s 2. q 2y I ^4fci+4fc 2 +4j. q 2^ 



y/(q 2 , g 4fc 1? ? 4fc 2j ? 4fc 1+ 4fc 2 +2j-2. gSy. y ( q 2 ]q 2) 



I n 



' 2 1 q \k. ■ I/,, - 1,. ;q ' q , 
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This determines C kl t ^ k n , but it takes some work to find the expression in terms of q- 
Hahn polynomials. First, take n± as the summation parameter in the sum and multiply 
both sides by to find that both sides are polynomials of degree n+j in s 2 . Apply 
[10, (III. 6)] to rewrite the 3</?2-series as a polynomial in s 2 and the (/-binomial theorem [10, 
(II.3)] to write {q 4k2 s 2 ;q 2 ) j as a polynomial in s 2 . Comparing next the coefficients on both 
sides gives an expression for the Clebsch-Gordan coefficients as a terminating 3 ^ 2 -series. 
To put it into the required form in terms of g-Hahn polynomials, we need to apply some 
transformations for 3 </? 2 -series, namely [10, (III. 13), (III. 11)]. The constant follows by a 
straightforward calculation. □ 

Combining Lemma 4.4 with the unitarity of the intertwining operator consisting of 
the Clebsch-Gordan coefficients results in the orthogonality relations for the g-Hahn and 
dual g-Hahn polynomials, cf. [10, §7.2]. 

We now have all ingredients to rewrite (4.11). Simplifying proves the following theo- 
rem. 

Theorem 4.5. With the notation (4.7) and (4.6) for the Al-Salam and Chihara polyno- 
mials and Askey- Wilson polynomials and x\ = fx(wi), = n{w2), n,j G Z +; fci, > 
we have 

n+j 



i(n-Z) 



n + j 
I 



QM- 2l ^^- 2 ^ Ak2 -^ri + 3 -q Z ) 



<i 2 



1=0 

x si(x 1 ;q 2kl w 2 ,q 2kl /w2\q 2 )s n+j -i(x 2 ;q 2k2 s,q 2k2 /s\q 2 ) = 
s n { Xl ;q 2k ^ 2k2+2 h, q 2k ^+ 2k2+2 i/s\q 2 ) Pj (x 2 ; q 2kl Wl , q 2k ' /w u q 2k2 s, q 2k2 /s\q 2 ). 

Remark 4-6. (i) Theorem 4.5 is a connection coefficient formula for orthogonal polynomials 
in two variables, orthogonal for the same measure, where the connection coefficients are 
given by the g-Hahn polynomials. The dual connection coefficient problem follows from the 
orthogonality for the Clebsch-Gordan coefficients or, equivalently, from the orthogonality 
relations for the dual g-Hahn polynomials. 

(ii) The case j = gives a simple convolution property for the Al-Salam and Chi- 
hara polynomials, since the g-Hahn and the Askey- Wilson polynomial reduce to 1. This 
was the motivation for Al-Salam and Chihara [2] to introduce the Al-Salam and Chihara 
polynomials as the most general set of orthogonal polynomials still satisfying a convolution 
property, see also Al-Salam [1, §8]. The case n = is also of interest, since then the g-Hahn 
polynomial can be evaluated and the Al-Salam and Chihara polynomial on the right hand 
side reduces to 1. In both cases we have a free parameter in the sum. 

(iii) Formally, in the representation space £ 2 (Z + )®£ 2 (Z + ) we have two bases of (gener- 
alised) eigenvectors for the action of Y S A, namely v k (x) and w(x\; X2). They are connected 
by Clebsch-Gordan coefficients, which are now expressible as Askey- Wilson polynomials; 

w{ } = y pM^q 2k ^ u q 2k ^/x u q 2k2 s,q 2k2 /s\q 2 ) 

The dual Clebsch-Gordan coefficient relation follows by integrating against the appropriate 
orthogonality measure for the Askey- Wilson polynomials. 



CONVOLUTIONS FOR ORTHOGONAL POLYNOMIALS 



21 



4.2. Racah coefficients and orthogonal polynomials. In the tensor product of three 
positive discrete series representations ttj £i <g> 7tk 2 <E> of U q (su(l, 1)) we have the same 
orthogonal bases as in §3.2 and we use the same notation as in (3.9)- (3.12). Similarly, we 
now have an intertwiner for the U q (su(l : l))-action in terms of (/-Racah coefficients; 



Again the constraints (3.14) hold. 

Proposition 4.7. Define 0: £ 2 (Z+) ® £ 2 (%+) ® £ 2 (%+) -> L 2 (R 3 , dm(x u x 2 , x 3 )) by 
^{e k n \®e k n \®e k n l)(x u x 2 ,x 3 ) 



= 5 ni (xi; q 2kl w 2 , q 2kl /w 2 \q 2 ) S n2 (x 2 ; q 2k2 w 3l q 2k2 /w 3 \q 2 ) S ns (x 3 ; q 2k H, q 2k3 /s\q 2 ), 



with the measure dm(xi, x 2 , x 3 ) given by 

dm(x l]q 2k ^w 2 , q 2kl /w 2 \q 2 ) dm(x 2 ; q 2k2 w 3 , q 2k2 /w 3 \q 2 ) dm(x 3 ; q 2k3 s, q 2kz /s\q 2 ), 

where Xi = /J,(wi). Then is a unitary map intertwining it / Cl ®7Tfc 2 ®7rfc 3 ((l <g> A)(A(Y S A)) 
with2(M Xl - f x{s))/(q- 1 -q). 

Proof. Observe that (1 <g> A)(A(Y S A)) = A 2 <g> A(Y S A) + Y S A <g> A(l). The proof now 
proceeds as the proof of Proposition 4.2. □ 

Proposition 4.8. (i) The following equality holds with Xi = fx(wi); 



XPj(x 3 ; q 2kl2 w 1 ,q 2kl2 /w 1 , q 2k3 s, q 2kz / s\q 2 )p J12 (x 2 ; q 2kl w u q 2kl /w u q 2k2 w 3 , q 2k2 /w 3 \q 2 ). 



( e (fci(fc 3 fc3)fc 33 )fc) ( Xl ? X2? X3 ) = C ., (A . 1? k 23 )C j23 (k 2 , k 3 ) S n ( Xl ; q 2k s, q 2k /s\q 2 ) 
xp r (x 2 ; q 2k 'w u q 2k '/w u q 2k23 s, q 2k23 /s\q 2 )p j23 (x 3 ; q 2k2 w 2 , q 2k2 /w 2 , q 2k3 s, q 2k23 /s\q 2 ). 



The constant Cj(k\, k 2 ) is defined in Proposition 4-3. 

The proof of Proposition 4.8 is slightly more complicated than the proof of its coun- 
terpart Proposition 3.12 due to the fact that we do not have a nice factorisation for © as 
in Remark 3.11. This is a consequence of the non-cocommutativity of the comultiplication 



Note that the occurrence of S n (xi; q 2k s, q 2k / s\q 2 ) on the right hand side corresponds 
to the intertwining property of the Racah coefficients as in the proof of the first statement 
of Proposition 4.3. 



(4.12) 




(ii) The following equality holds with Xi = fji(wi); 



for C/ 5 (su(l,l)). 
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Proof. The proof of (i) and (ii) is similar. To prove (i) we use (3.10) (for the U q (su(l, In- 
setting) and Proposition 4.7 to find 

Q(e« k ^ k ^ k )(x 1 ,x 2 ,x 3 ) = Yl ^:;«^«3fe;9 a3 M 2t V^ 2 ) 

ni2+n 3 =n+j 

X E C knl%\ 2 2 S »! ^2, ^AH^) ^n 2 (x 2 ; ^3, tf 1 * 

ni+n 2 =n 12 +ji2 

=C jl2 (k u k 2 )p jl2 (x 2 ; q 2k 'w u q 2k ^/w u q 2k2 w 3 , q 2k * /w 3 \q 2 ) 

X E C n\V k nfn Sn 3 (x 3 ; q 2k *s, q 2k * / s\q 2 ) S ni2 ( Xl ; q 2k ™w 3 , q 2k ^/w 3 , \q 2 ) 
n 12 +n 3 =n+j 

by (4.11). The last sum can be evaluated by another application of (4.11) leading to the 
result. □ 

The n-dependence in the right hand sides of Proposition 4.8 is the same, so we obtain 
the Wigner-Eckhart theorem for the U q (su(l, l))-setting by applying O to (4.12). So we 
can restrict to n = in (4.12) before applying without loss of generality, and we obtain 
the following polynomial identity in x 2 and x 3 with w\ as a parameter; 

(4.13) Cj(k 12 , h)C J12 (h, k 2 ) Pj (x 3 ; q 2kl2 Wl , q 2k ^/ Wl , q 2k *s, q 2k3 /s\q 2 ) 
x p jl2 (x 2 ; q 2k ^w u q 2k ^/ Wl , q 2k *w 3 , q 2k */w 3 \q 2 ) = ^ <Vfi£, ia C i'^ k 23 )C hs (k 2 , k 3 ) 

J23=0 

x Pj ,(x 2 ; q 2kl w u q 2kl /w u q 2k23 s, q 2k23 /s\q 2 ) p hs (x 3 ; q 2k2 w 2 , q 2k2 /w 2 , q 2k H, q 2k23 /s\q 2 ). 

Again we can use (4.13) in two ways. Firstly, we specialise to a suitable formula from 
which the Racah coefficients can be determined explicitly. Secondly, with the explicit 
expression for the Racah coefficients we derive a convolution identity for the Askey- Wilson 
polynomials. 

Proposition 4.9. The Racah coefficients of (4. .12) are given by 

, , , f f Ak 1 +4k 2 +2j 12 -2 r Ak 2 +4k 3 +2j 23 -2 n ~2j 12 n ~2j 23 \ 

jjkiMMi _ f> I y >y >y ■n 2 -n 2 \ 

u k 3 ,k,k 2 3 U4W 1 q 4k 2 ^ q 4k 1 +4k 2 +4k 3 +2j+2j 12 -2^ q -2j-2j 12 I) 

wit/i t/ie constant C given by 

( n 2 4fci „4fc 23 „4fci+4fc 2 3 +20-+ji2-j23)-2. „2\l/2 / 2 4fc 2 „4fc 3 4fc 2 +4fc 3 +2j 23 -2. 2\l/2 

^ ^ ^ ^ >y )j+ji2-j 23 y L i >y >y ^23 

(Q 2 , g4fcia, g4fcs } g4fci 3 +4fc 3 +2j-2. g2)l/2^ 2j ^ ? ? 4fc 2 ? ? 4fc 1+ 4fc 2 +2j 12 -2. g 2)V2 

(? 4fc3 ;? 2 ) J (? 4fc2 ;? 2 ) J12 (? 4fcl+4fc2+4fc3+2j+2jl2 - 2 ;? 2 ) J23 



x q2k 2 (j-j 23 ) 



J +J12 
323 



72 (Q 4fc 3,^ 2 +4/c3+2 j2 3-2. ? 2 )j23(g 4 fc2+ 4 fc 3+4 j2 3 ;? 2 )j+ji2 _ j23 - 



Proo/. Take wi = s = 1 and x 2 = ^(q 2kl ) in (4.13) to find 

Pj(x 3 ; q 2k ^\q 2k ^\q 2k ^q 2k3 \q 2 )(q 2k ^ 2k2 w 3j q 2k ^ 2k * /w 3 ; q 2 ) jl2 

323, 
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for C\ an explicit constant depending upon fci, hi, k 3 , j'12,.7'23 an d J, since two Askey- 
Wilson polynomials can be evaluated for this choice. So the Racah coefficients occur as 
the coefficients when developing the polynomial of degree j + ji2 on the left hand side into 
Askey- Wilson polynomials. Hence the Racah coefficients can be obtained from 



U2U k 3 ,k,k 2S 



Pj23 ( x 3; q 2kl+2k2 , q 2k 2- 2k i ? q 2k 3 ^ q 2k z \q 2 ^p- ( X3 ; q 2kl2 , q 2kl2 , q 2kz , q 2ks \q 2 ) 

x (q 2k ^+ 2k2 w 3 , q 2k ' +2k2 /w 3 ; q 2 ) J12 dm(x 3 ; q 2kl+2k2 , q 2k2 ~ 2kl , (p 3 , q 2k " |(/ 2 ), 

for some known constant Ci- Now observe that 

(q 2k ^ 2k2 w 3 , q 2k ^ 2k2 /w 3 - q 2 ) J12 dm(x 3 ; q 2k ^ 2k2 , q 2k2 ~ 2k ^ , q 2k * , q 2k * \q 2 ) 

= C 3 dm{x 3 - ,2*12^2*2-2*! ; q 2k^ q 2k 3 ^ 

for some constant known C 3 by (3.14) and (4.8). Thus the Racah coefficients can be 
obtained by integration; 

(4-14) C 4 u k k ^ = Jp J23 (x 3 -q 2k ^ 2k2 iq 2k2 - 2k \q 2k ^q 2kz \q 2 ) 

x Pj (x 3 ; q 2kl2 , q 2kl2 , q 2k3 , q 2k * \q 2 ) dm(w 3 ; q 2kl2 , q 2k2 ~ 2k ^ , </ 2fc 3 ? q 2k s \ q 2 )^ 

with C 4 explicitly known. Observe that three out of four of the parameters of each of 
the Askey- Wilson polynomials in (4.14) coincide with the parameters of the Askey- Wilson 
measure in (4.14). Use the connection coefficient formula for Askey- Wilson polynomials 
with one different parameter, cf. Askey and Wilson [3, (6.4-5)] or see [10, (7.6.8-9) with 
the right hand side of (7.6.9) multiplied by (q;q) n ], twice to rewrite the Askey- Wilson 
polynomials in terms of Askey- Wilson polynomials with the same parameters as the Askey- 
Wilson measure in (4.14). By orthogonality the integration is then easily performed and 
we are left with a single sum, which can be written as a very-well-poised 8</ : '7" ser i es - This 
can be transformed to a balanced 4</?3-series by Watson's transformation [10, (III. 17)], 
and another application of Sears's transformation [10, (III. 15)] gives the form as in the 
statement of the proposition. The constant follows from bookkeeping. □ 

Recall the (/-Racah polynomials, see [10, §7.2], [15], 

/ q~ n , a(3q n+1 , q~ x , -f5q x+1 
(4.15) Rn(i/(x);a,l3,^,6;q) = 4^3 [ ' a \ ;q,q 

with v(x) = q~ x + ^5q x+1 : one of the lower parameters equals q~ N , N e Z+, and < n < 
N. The 4^3-series in Proposition 4.9 can be written in terms of a (/-Racah polynomial. 

We can now rewrite (4.13) to arrive at the key result of this paper. For convenience 
we replace q 2 by (/, (a, 6, c) = (q kl , q k2 , (/ fcs ), and we relabel w\, j'23 and j\2 by t, I and n, 
and finally replace X2, x 3 by x±, X2- We obtain the following (/-analogue of Theorem 3.13 
and Corollary 3.15. 
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Theorem 4.10. With x\ = fi>(wi), x 2 = n(w 2 ), n,j G Z + , we have the convolution 
identity for the Askey-Wilson polynomials 

n+j 



5>- 



1=0 



j +n 
I 



(b 2 ;q) n (a 2 b 2 c 2 qi+ n - 1 ;q) l (c 2 ;q) 



q (c 2 ,b 2 c 2 q l 1 ;q)i(b 2 c 2 q 2l ;q) j+n -i 
x R^n); b 2 /q, c 2 /q, q~^~ n ~ 1 , oV?^" 1 ;?) 
x pj +n _i(xx; at, a/t, bcq l s, bcq 1 / 's\q) pi(x 2 ; bw 1 , b/w 1 , cs, c/s\q) 
= p n (xi,at, a/t, bw 2 , b/w 2 \q) Pj(x 2 ; abq n t, abq n /t, cs, c/s\q), 
with the notation of (4.6), (4.15). 

Remark 4-11- (i) Theorem 4.5 can be obtained as a special case of Theorem 4.10 by letting 
0. 

(ii) Theorem 4.10 leads to a kind of generating function for the g-Racah and g-Hahn 
polynomials. Choosing w\ = at and w 2 = cs in Theorem 4.10 reduces all four Askey- 
Wilson polynomials to a single term. The remaining free parameters s and t appear only 
in the combination s/t. Replacing bcs/(at) by u, and (a 2 ,b 2 ,c 2 ) by (a, (3, 7) gives: 

(a;q) j+n -i(/3;q) n (a/3jq J+n - 1 ;q)i 



E 



1=0 



j +n 
I 



-v? l ((3-fq l /u;q) j+n -i(u;q)i 



q (a;q)n(Piq l 1 ; q)l(Piq 2L , q)j+n-i 

x R l (v(n);l3/qn/q,q- i - n -\<xl3q n+j - 1 -,q) = (aq n u;q) 3 {P /u;q) n . 
For 7 = we obtain a similar identity for g-Hahn polynomials: 

n+j 



E 

1=0 



j + n 
I 



(<x;q)j+n-i(P; q)r 

(ar,q)n 



Qn(q l ;(3/q,a/q,n+j;q)u J L (u;q)i = (aq n u;q)j((3/u;q) n . 



.3-1 



(iii) Theorem 4.10 gives the connection coefficients for two sets of orthogonal polyno- 
mials with respect to the absolutely continuous measure 



±2 „„±2 



wr-. 



! ;?)c 



1 



dwi dw 2 



(taw 1 ,aw 1 /t,csw 2 ,cw 2 /s;q) OQ (bw 1 w 2 ;<?)oo w i w i 

on the torus T 2 for < \a\ < \t\, < |c| < \s\ and \b\ < 1. Here all possible signs for 
± have to be used. (Otherwise discrete masses at points and lines have to be added.) This 
weight function is invariant under simultaneously interchanging w\ with w 2 , t with s and 
a and c. This transforms the orthogonal polynomials on the right hand side of Theorem 
4.10 to the ones occurring in the left hand side. 

In particular, note that for s = t, a = c, the weight function is invariant under the Weyl 
group for B 2 , i.e. the group generated by {w\,w 2 ) 1— > (w 2 ,w\) and ( w\,w 2 ) 1— > (wi, w^ 1 ). 
The corresponding Weyl group invariant orthogonal polynomials are 

p n (fi(wi); at, a/t, bw 2 , b/ w 2 \q) pj(/i(w 2 ); abq n t, abq n /t, at, a/t\q) 

+ p n (fi(w 2 ); at, a/t, bw\, b/w\\q) pj{p,{w\); abq n t, abq n /t, at, a/t\q) 

for n > j > 0. These orthogonal polynomials do not seem directly related to the 
Koornwinder-Macdonald orthogonal polynomials associated with root system BC 2 , see 
[20], although the structure of the orthogonality measure is similar. 
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5. Linearisation coefficients for Askey- Wilson polynomials 

In the results of the previous section using U q (su(l, 1)), especially Theorems 4.5 and 
4.10, we can use analytic continuation with respect to the parameters involved to find 
similar identities but with the Al-Salam and Chihara polynomials and the Askey-Wilson 
polynomials replaced by the dual g-Krawtchouk polynomials and the g-Racah polynomials. 
These identities can be obtained by the same procedure using U q (su(2)) and its represen- 
tation theory instead of using U q (su(l, 1)). In particular we can now give an interpretation 
for g-Racah polynomials as Clebsch-Gordan coefficients for U q (su(2)). In this case we also 
have some knowledge on the structure of the dual Hopf *-algebra and this can be used to 
obtain a linearisation formula for a two-parameter family of Askey-Wilson polynomials. 
This is an application of the results of the previous section. 

We first recall U q ($u(2)) and its representation theory, see e.g. [6], [17], [21]. The Hopf 
algebra structure on U q (su(2)) is the same as the Hopf algebra structure on U q (sl(2 : C)), 
cf. (4.1), (4.2). The *-operator making U q (su(2)) into a Hopf *-algebra is given by 

A* = A, B* = C, C* =B, D* = D. 

There is precisely one irreducible unitary L r g (su(2))-module Wn of each dimension 
N + 1 with highest weight vector v +: i.e. Av + = q N ^ 2 v +1 B v + = 0. The corresponding 
representation is denoted by t N . With respect to the standard ortho normal basis , 
< n < N, the action of the generators is given by t N (A) e% = q n ~ N / 2 e% and 



t N (B) el = q -Y^r\J{l - Q 2n+2 )(l ~ Q 2N ~ 2n ) e% +1 , 
t N (C) = q - Y —^\J( 1 ~ ^)(! " Q 2N ~ 2n+2 ) eJLx, 



with the convention = = ej^ +1 . So e 1 ^ is the highest weight vector. The represen- 
tation t N , considered as a representation of U q (sl(2, C)) can be obtained from the discrete 
series representation tv^ of (4.3) by formally replacing k by —N/2. 

The Clebsch-Gordan decomposition holds; as unitary L r (? (su(2))-modules 

min(AT 1 ,iV 2 ) 

w Nl ®w N2 = W N 

3=0 

The matrix coefficients of the intertwining operator give the Clebsch-Gordan coefficients; 



(5-1) ^=^CWe>e5 

m,n 2 



Of course, the Clebsch-Gordan coefficient is zero if iV ^ N\ + N2 — 2j for < j < 
min(A"i, N 2 ). By considering the action of A on both sides we see that the Clebsch-Gordan 
coefficient is zero unless n\ + n 2 = n + j, so that the sum is actually a single sum. The 
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Clebsch-Gordan coefficients are normalised by (e^, e^ 1 £g> e^ 2 ) > if N = Ni + N 2 — 2j, 
< j < min(AT 1 , N 2 ). 

We are particularly interested in the element 

X p = q 1 ' 2 B + q~ x l 2 C - pl/2 q Z P q -i\ A ~ D ) e U q (sl(2, C)), p > 0. 

Then X P A is self-adjoint and A(X P A) = A 2 ®X P A+X P A®1. Koornwinder [21] has shown 
that in each module Wn the action of X P A is completely diagonalisable. To formulate this 
result we introduce the orthonormal dual g-Krawtchouk polynomials; for a > 0, 



r n (q~ x ~ q x ~ N /a; a, AT; q) = (-l) n a n/ V (n ~ 1)/4 



N 
n 



1/2 



R n {q- x -q x - N /a;a,N;q), 

Q 



o / ,-/v / v fq n ,q x ^-q x N / a 

R n {q -q /a; a, N; q) = 3 tp 2 ( _ N ;q,q 



for N G Z_|_ and < x, n < N. The dual g-Krawtchouk polynomials are defined by 

' ) q 

q~ N ,0 

The corresponding three-term recurrence relation is 

(q~ x - q x ~ N /a) r n = A n r n+1 + q n ~ N (l - a" 1 ) r n + A n _i r„_i, 

(5.2) 



A n = a -V2 ? -iV+n/2 v / (l _ q n+l ){1 _ q N-n^ 

for < x, n < N and r n = r n (q~ x — q x ~ N /a; a, N; q). 

Proposition 5.1. ([21]) There exists an orthogonal basis ^ = 4>J (p), < f < N, of 
Wn of eigenvectors of t N (X p A) for the eigenvalue 

,7V, N p l/2 q N-2f _ p -l/2 q 2f-N +p -l/2 _ p l/2 
X f (P) = =1 • 

q q 

Moreover, (f>f (p) = Yln=o r n(q~ 2f - q 2f ~ 2N /p;p, N;q 2 ) e% . 

Proposition 5.1 is the analogue of Proposition 4.1, and we could have formulated it in 
a similar fashion using the finite discrete orthogonality measure for the dual g-Krawtchouk 
polynomials. Actually, replacing e l6> , k in S n (cos9; q 2k s, q 2k / s\q 2 ) by q~ 2 f +N j 's, —N/2 and 
next taking s 2 = — p~ x gives r n (q~ 2 ^ — q 2 f~ 2N /p;p 7 N;q 2 ). The analogue of Proposition 
4.2 is the following. 

Proposition 5.2. For < fx < Ni, < f 2 < N 2 define in Wn x <E> Wn 2 the vector 

& = £ r ni (q- 2h - q 2h - 2N '- 2N2+4h /p;pq 2N2 - 4h ,N i; q 2 ) e% 
m=0 
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then t Nl <g) t N2 (A(X p A)) 0^j^ 2 = >$£g a (p) ■ Moreover, < f x < N lt 

< f 2 < AT 2 , constitutes an orthogonal basis o/Wjvj <E> Wjv 2 o/ eigenvectors of A(X P A) . 

Proof. From A(X P A) = A 2 <g> X P A + -X" P A <g> 1 it follows that there is an eigenvector of the 
form ^^=0 Pm e nt ® ^ 2 by solving a three-term recurrence relation for the p ni . Now 
(5.2) can be used to solve this. 

There are {N x + 1)(N 2 + 1) eigenvectors in W Nl <g> W N2 and {^j 2 2 Agl^) equals 
zero if f 2 ^ gi by Proposition 5.1 and it also equals zero if fi + f 2 ^ gi + gi by the 
self-adjointness of X P A. □ 

The result of Proposition 5.2 can be obtained from Proposition 4.2 by substituting 
fci, k 2 by —Ni/2, -N 2 /2 and w 2 by q N ^+ N *- 2 h ~ 2 h/ s , q N ^~ 2 ^ / s and s 2 by -p _1 . 

Proposition 5.3. For JV = JVi + iV 2 - 2j, < j < min(AT 1 , N 2 ) we have 
if < fi + f 2 — j < N, and (^/^/^ e n) = otherwise. If non-zero, then 



A^2 

J 



1/2 j/2 j{2N 1 +N 2 ) _3j0'-1)/2 

g 9 9 C -1 2/ 1 +2/ 2 -2Af 1 -2iV 2 2^ . 

g2 y/(q^i jq 2N 1+ 2N 2 -2j+2. q -2^ F 

/ -2j 2j-2-2N x -2N 2 -2/ 2 _„-l„2/ 2 -2JV 2 \ 
/ -2/i-2/ 2 2\ / 9 >9 ,9 , P 5 2 2 \ 

x ^ ,g j j4 ^3 ^ q _ 2N2 ^ q _ 2fl _ 2f2 ^_ p _ lq2fl+2f2 _ 2Nl _ 2N2 ,q ,q y 

Note that the 4</?3-series is balanced, and can be written in terms of g-Racah polyno- 
mials (4.15). The 4</?3-series in Proposition 5.3 equals 

i2 . (? -2/ 2 _ p -l ? 2/ 2 -2iV 2;? -2iV 2 -2 j? -2Ar 1 -2 j? -2/ 1 -2/ 2 -2 j _ p -l ? 2/ 1+ 2/ 2 -2AT 2 . ?2) _ 

The proof of Proposition 5.3 is similar to the proof of Proposition 4.3. Proposition 5.3 can 
also be obtained from Proposition 4.3 using the substitutions as indicated earlier. It can 
also be obtained by using = Yl ^nin^o* e nl ® e n 2 : Propositions 5.2 and 5.1 and the 
explicit value for the Clebsch-Gordan coefficients for n = 0, N = Ni + N 2 — 2j, 



r N 1 ,N 2 ,N =( _^ n2n n 2 (N 2 - j -l) I (q 2ni+2 ; q 2 )n 2 (g 2Ar i~ 2n i ; Q~ 2 ) n2 / (g^g^jj 

ni '" 2 ' [ } Q V (<l 2 ^ 2 )n 2 (q 2N2 ;q- 2 )n 2 " V {q^ + N 2 -2 3+ 2. q 2 )3 - 

See e.g. [24, §14.3], but this simple case can also be derived as follows. Apply C to both 
sides of (5.1) to obtain a three-term recurrence for the Clebsch-Gordan coefficients, which 
reduces to a two-term recurrence for n = 0. This can be easily solved, with the initial 
condition following from the unitarity and the normalisation, see [23] for a similar deriva- 
tion. Then we have a sum involving the product of two dual g-Krawtchouk polynomials. 
Upon inserting the series representation we obtain a triple sum, and after interchanging 
summations we can use the (/-binomial theorem and the g-Chu-Vandermonde sum, see [10], 
to obtain a single 4 </? 3 -series. 
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Remark 5.4- With Proposition 5.3 at hand it is straightforward to calculate the U q ($u(2))- 
counterparts of Theorems 4.5 and 4.10. The Al-Salam and Chihara, respectively Askey- 
Wilson, polynomials have to be replaced by dual g-Krawtchouk, respectively g-Racah, 
polynomials. The result can also be obtained from Theorems 4.5 and 4.10 by substitution 
as indicated earlier, so we do not give them explicitly. These formulas give an alternative 
for the formulas of Groza and Kachurik [12]. 

In the representation space Wn 1 <E> Wn 2 we have two bases of eigenvectors for the 
action of X p A, namely 4>J and (p^fe 2 , an d the corresponding Clebsch-Gordan coefficients 
are given by Proposition 5.3, since, with N = Ni + N 2 — 2j, 

min(Ar l5 Ar 2 ) N 

&= E E<*XT.O«? 

j=0 n=0 
min(iVi,iV2) N 

= E ,e»)J2r n (<l 2j - 2h - 2h - q 2h+2h - 2 >- 2N /p;P,N ]q *)e» 

j=0 n=0 
min(N 1 ,N 2 ) 

V (fh Nl ' N2 P N \rh N 

j=o 

Here we use the convention that = for / > N or / < 0. So introducing the notation 

(s-3) Of = E ^& N (p) tf. 

f,j 

we see that the Clebsch-Gordan coefficients are zero unless fi + f 2 = / + j, and then 
C^fa*f N (p) = (0/^/^ 2 ) e o)- So, by Proposition 5.3 we have proved that the g-Racah 
polynomials occur as Clebsch-Gordan coefficients for U q (su(2)). 

Using (5.3) in a special case we can obtain the linearisation coefficients for the two 
parameter family of Askey- Wilson polynomials occurring as spherical functions on the 
quantum SU(2) group, cf. [21]. We consider odd-dimensional representations; Ni = 2li, 
N 2 = 2/ 2 , h,l 2 e Z+. Then the kernel of t 2ll (X p A) is one dimensional and spanned by 
(j) 2l *{p). Moreover, 0f Vf 2 ^) = ^i^ip) ®4> 2l2 {p)- Next we consider matrix elements as 
linear functionals on U q (su(2)) to find 

(5.4) = ( t2h ®t 2/2 (A(X))^f(p),^;f 2 (r)> 

= E ^ff ,2i b)^; 2 ^ 2 'W(^(^)^(p),^(r)), 

i=\h-h\ 

where r > is another parameter and A(A) = J2(x) -^(i) ® -^(2)- 
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The dual Hopf *-algebra A q (SU(2)) generated by the matrix elements of the repre- 
sentations t N , N E Z_|_, of £/" g (su(2)), is known in terms of generators and relations, cf. [6], 
[17], [21]. Koornwinder [21] has given an explicit expression for the element in A q (SU(2)) 
corresponding to the linear functional considered in (5.4); 

(?\X)<j>?{ P U?{r)) = {q2l %[ q2)i (XMP-,^,^,^,-^^)), 

where p G A q (SU(2)) is some fixed simple element, which is, up to an affine scaling, the 
linear functional X i— > {t 2 {X)(j)\ (p), <f>\ (r)), and the last (•,•) denotes the duality between 
U q (su(2)) and A q (SU(2)). Since A q (SU(2)) is the dual Hopf *-algebra, the left hand side 
of (5.4) corresponds to the multiplication of the two linear functionals. So (5.4) leads to 
the following identity in A q (SU(2)); 

h+l2 ( n 2h+2. „2\ /„2Z 2 +2. _2\ 

i=\h-h\ {q ' q )l 

with pi(-) = p,(.; ?v /f , Qy^", ^L, -q^p¥\q 2 ). 

The only information on A q (SU(2)) needed is the existence of a family of one- 
dimensional representations sending p to cos 9. Thus, applying the one-dimensional rep- 
resentations of A q (SU(2)) and using Proposition 5.3 proves the following linearisation 
coefficient formula. 

Theorem 5.5. Let pi(x) = pi(x; Qa/t", 1\J^^ ~Jpr J ~Q.\fw\ c L 2 )j P,r > 0, be defined in terms 
of Askey- Wilson polynomials (4.6). Then the coefficients in the linearisation formula 

2 min(7i,Z 2 ) 

p h (x) p h (x) = E C 3 Pli+h-j( X ) 
3=0 

are given by a product of two balanced terminating ^ip^-series; 

^2/i+2/ 2 . q~2y l _ ? 4/i+4/ 2 -4j+2 

X ( q 2h+2h+2. q 2 ) h+h _ J 1 - ^i+4Z 2 -2 J+ 2 

f a~ 2 3 n~ 212 „2j-2-4Ji-4«2 _„-l„-2« 2 \ 
,72/ -1 -2L-21? 2\ / */ > */ '*/ ■> P H 221 

xr (" r ^ 2 ;y) j4 ^ 3 ^ ? _ 4Z2 ^_ 2Zl _ 2Z25 _ r _ 1? _ 2Zl _ 2Z2 

Remark 5.6. (i) In particular, for p = r the linearisation coefficients are positive. This can 
already be observed without the explicit knowledge of the linearisation coefficients, see [17, 
§8.3], [19, §7]. 
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(ii) For p = r = 1 the Askey- Wilson polynomials pi(x; q, q, —q, — q\q 2 ) are the continu- 
ous g-Legendre polynomials Ci(x; g 2 |g 4 ), see [3, §4]. This is a special case of the continuous 
g-ultraspherical polynomials introduced by Rogers at the end of the 18th century. Rogers 
calculated the linearisation coefficients for the continuous g-ultraspherical polynomials, see 
e.g. [3, §4], [10, §8.5], and we can go from Theorem 5.5 to the special case of Rogers's 
result by using Andrew's summation formula, see [10, (11.17)]. 
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